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Abstract 

In this talk I compare the knowledge on the critical end point of the QCD phase diagram grasped 
from lattice calculations, with that obtained from Nambu-Jona-Lasinio (NJL) model computations. 
The original publication is available at ' |http : / / www. sif . it/SIF/en/p ort al / j ournals] ' . 



The major knowledge on the QCD phase transitions at zero baryon density comes from 
first principle calculations made on supercomputers, namely from the lattice. When simula- 
tions are run with physical quark masses, it is well known that lattice predicts the restoration 
of chiral symmetry, which is spontaneously broken by the quark condensate in the vacuum, 
at a finite value of the temperature 170 MeV < T < 200 MeV. The chiral restoration in 
the vacuum is actually a smooth crossover, the reason being that finite values of the quark 
masses break explicitly chiral symmetry, hence there is not a true phase transition. For sim- 
plicity, from now on I will call the chiral crossover, as well as the true phase transition, the 
chiral restoration. In correspondence of the chiral restoration, lattice shows that a decon- 
finement transition occurs. This has suggested that chiral restoration and deconfinement of 
color are two intimately connected transitions of QCD, (see Ref. [l[ and references therein). 

Lattice investigations at finite baryon chemical potential, /i, suffer the (in)famous sign 
problem in three color QCD. To this end, several approximated methods have been used to 
circumvent it. By means of one of these methods, namely the two parameter reweighting, 
it has been predicted, some time ago {2J], that the chiral crossover becomes a first order 
transition at a certain value of /i. The couple (fiE, Te) in the (/i,T) plane at which this 
occurs is called the Critical End Point (CEP) of the QCD phase diagram. The numerical 
simulations of Ref. [J predict fi E ~ 350 MeV and T E « 160 MeV. 

An interesting alternative to the reweighting analysis of the QCD phase diagram, with 
particular reference to the existence of a CEP, has been performed p] (see also references 
therein). The reasoning on which the investigations of js} lies is very simple to understand: 
at /i = 0, it is known, from lattice studies, that the chiral transition is a true first order 
transition, if quarks are taken in the chiral limit. Moreover, as the quark masses are increased 
above a critical value, the transition becomes a crossover. It happens that at the physical 
point, defined as the couple of values for the up- and strange-quark mass, (m u ,m s ), which 
gives the physical spectrum of mesons, the transition is a crossover. Hence, there exists a 
critical line in the (m u , m s ) plane which is the border between an inner region, in which the 
chiral transition is of first order, and an outer region, in which the transition is a crossover. 
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the Taylor expansion, the critical line is expressed as 
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The coefficient c\ governs the behavior of the critical line at small values of \i. Nowadays, 
the coefficients Ck have been determined up to the 8 th order. Surprisingly enough, the results 
of Ref. 3| are that the critical line moves towards lower values (hence to less realistic) of 
the quark masses, as \x is increased. This means that at finite (but small, see below) values 
of \x the crossover remains crossover, if quarks are taken in the chiral limit. The analysis 
performed in Ref. [s| should be reliable, by author's admission, up to fi ~ 500 MeV. As a 
consequence, their results are consistent with the scenery in which a CEP, if it there exists, 
is located at values of ft larger than that predicted in |2j. The discrepancy is probably due 
to the fact that the reweighting method suffers of large systematic errors at large //. 

It is of a certain interest to compare this scenario with that of some model calculation. 
Among the various models, the NJL model (or its improved version, the Polyakov-Nambu- 
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Jona-Lasinio (PNJL) one) is a very popular one (for review see 
grangian shares the same global symmetries of the QCD Lagrangian. Since we can describe 
the numerous (expected) phases of the QCD phase diagram in terms of broken/restored 
global symmetries, the hope is that the NJL calculations grasps, for the property specified 
above, at least the main characters of the QCD phase diagram in the \x-T plane. Moreover, 
determining the ground state of the model at any temperature and/or chemical potential 
is a very easy task, which requires only some numerics. On the other hand, first principle 
calculations are not feasible at finite \x both numerically (for the infamous sign problem of 
three color QCD) and analytically (for weak coupling approximation might break down in 
the range of temperature/chemical potential relevant for heavy ion collisions as well as for 
compact star phenomenology). Therefore, the NJL model might be helpful in depicting the 
main aspects of the QCD phase diagram. 

The NJL (or PNJL) phase diagram has been discussed in several papers. Here I refer 
to (5)]. First of all, I need to specify the model Lagrangian density, 

C = Y,^f( W ^- m f + ^o)^f + A + £ 6 , (2) 
/ 

where the sum is over the three flavors / (= 1,2,3 for u,d,s). In the above equation the 
background gauge field = g5^A a ^T a is coupled to quarks via the covariant derivative 



Dfj, = — iAfj_ and will be specified later; m/ is the current mass (we assume m u = 
ma)- The quark chemical potential is denoted by fi. The NJL four-fermion and six-fermion 
interaction Lagrangians are as follows 4j: 

8 

C 4 = GJ2 [(^) 2 + (#7 5 A^) 2 ] , (3) 

a=0 

£ 6 = -K[deti>f{l + l5 )ip f + deti) f (l- l5 )i) f ] , (4) 



where A a are the Gell-Mann matrices in flavor space (Ao = a/2/3 1/) and the determinant 
is in flavor space as well. The parameters are 

m ud = 5.5MeV, m s = 140.7 Me V, GA 2 = 1.835, KA 5 = 12.36, A = 602.3 MeV. 



From these parameters one gets m n ~ 135 MeV, mx — 498 MeV, ~ 958 MeV, m 
515 MeV and / w ~ 92 MeV. 

Once the Lagrangian is specified, the thermodynamic potential at temperature T is ob- 
tained after integration over the fermion fields in the partition function: 

n = w[r,$,$] + n,[M / ,$,$] , (5) 

where f2 g denotes the free quark contribution, as well as the interaction term of quarks 
with the Polyakov loop (see [s| for more details). In the thermodynamical potential, the 
term U(T, $, $) is the novelty that improves the NJL model and promotes it to the PNJL 
model [6J. It describes the dynamics of the traced Polyakov loop in absence of dynamical 
quarks. The potential U cannot be determined by first principles: one has to chose a conve- 
nient form for it, by trying to reproduce lattice data on thermodynamical quantities of the 
pure glue theory. Different analytical forms of U lead to different quantitative predictions, 
even if the qualitative picture is quite not sensible of the form chosen. In this talk I focus 
on a model calculation based on the following potential: 

U<yT ^ $) = + b(T) ln[l - 6$$ + 4($ 3 + $ 3 ) - (<S$) 2 ] (6) 

where the analytical form of the coefficients has been determined in Ref. 0]. 

In the mean field approximation, which is formally equivalent to determine only the clas- 
sical contribution to the partition function, one can get quark condensates cry and Polyakov 
loop for any value of the parameters \i and T simply by looking at the global minima of Q. 
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FIG. 1: Sketch of the phase diagram of the PNJL model with 2+1 massive flavors. Here u q denotes 
the quark chemical potential, u q = 3u, where u corresponds to the baryon chemical potential. For 
simplicity, I have drawn only the chiral crossover line. The dashed line denotes the chiral crossover, 
the solid line corresponds to a first order transition. The region denoted symbolically by XSB 
denotes the zone of the phase diagram with quark condensate different from zero. In the region 
X ~ 0, on the other hand, one has (uu) m but (ss ^ 0). Based on 

Depending on the values of cr/ and <3>, one can characterize the symmetry breaking pattern 
of the theory in any point of the plane fi-T, hence one can build a phase diagram. The 
phase diagram of the model is sketched in Fig. [TJ For simplicity, I have drawn only the 
chiral crossover line. The dashed line denotes the chiral crossover, the solid line corresponds 
to a first order transition. The region denoted symbolically by xsb denotes the zone of the 
phase diagram with quark condensate different from zero. In the region x ~ 0, on the other 
hand, one has (uu) ~ but (ss) ^ 0. It is interesting to notice that, with the parameters at 
hand that reproduce the vacuum spectra of the pseudoscalar mesons, the CEP is located at 
quite large values of the quark chemical potential, which is one third of the baryon chemical 
potential, thus at values of /i larger than the 350 MeV quoted above. The introduction of a 
vector interaction can shift /xcep to higher values, depending on its magnitude at finite den- 
sity jsj]. It can even disappear at all, if the vector interaction is repulsive enough. Hence, we 
can conclude that the PNJL model scenario is in agreement with the newest lattice findings 
on the absence of a CEP at small values of the baryon chemical potential. 



5 



Acknowledgments 



I acknowledge H. Abuki, M. Ciminale, P. Colangelo, R. Gatto, M. Mannarelli and S. 
Nicotri for valuable collaboration and interesting discussions on the topics discussed here. 
The original publication is available at ' |http: / / www.sif.it / SIF /en/ portal/journals| ' . This 
talk is dedicated to the memory of my mentor, prof. Giuseppe Nardulli. 



[1] A. Bazavov et al, arXiv:0903.4 379 [hep-lat]. 

[2] Z. Fodor and S. D. Katz, JHEP 0404, 050 (2004) [arXiv:hep-lat /0402006] . 

[3] P. de Forcrand and O. Philipsen, PoS LATTICE2008, 208 (2008) [arXiv:Q811.3858l [hep-lat]]. 

[4] S. P. Klevansky, Rev. Mod. Phys. 64, 649 (1992); T. Hatsuda and T. Kunihiro, Phys. 

Rept. 247, 221 (1994) l; arXiv:hep-ph/9401310l ; M. Buballa, Phys. Rept. 407, 205 (2005) 

[arXiv:hep-ph/0402234] . 

[5] M. Ciminale, G. Nardulli, M. Ruggieri and R. Gatto, Phys. Lett. B 657, 64 (2007) 

|arXiv:0706.42T5l [hep-ph]] . 
[6] K. Fukushima, Phys. Lett. B 591, 277 (2004) [arXiv:hep-ph/0310121| . 

[7] S. Roessner, C. Ratti and W. Weise, Phys. Rev. D 75, 034007 (2007) |arXiv:hep-ph/0609281| . 
[8] K. Fukushima, Phys. Rev. D 77, 114028 (2008) [Erratum-ibid. D 78, 039902 (2008)] 
|arXiv:0803.33T8l [hep-ph]]. 



6 



